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Abstract 

We consider cascade models of turbulence which are obtained by restricting 
the Navier-Stokes equation to local interactions. By combining the results of 
the method of extended self-similarity and a novel subgrid model, we investi- 
gate the inertial range fluctuations of the cascade. Significant corrections to 
the classical scaling exponents are found. The dynamics of our local Navier- 
Stokes models is described accurately by a simple set of Langevin equations 
proposed earlier as a model of turbulence [Phys. Rev. E 50, 285 (1994)]. This 
allows for a prediction of the intermittency exponents without adjustable pa- 
rameters. Excellent agreement with numerical simulations is found. 
PACS-numbers: 47.27.Jv, 02.50.Ey, 47.27.Eq, 47.11. +j 
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I. INTRODUCTION 



Much of our intuitive understanding of turbulence is based on the concept of interactions 
which are local in k-space. Physically, it is based on the notion that most of the distortion 
of a turbulence element or eddy can only come from eddies of comparable size. Turbulent 
features which are much larger only uniformly translate smaller eddies, which does not 
contribute to the energy transfer. This immediately leads to the idea of a chain of turbulence 
elements, through which energy is transported to the energy dissipating scales. Accepting 
such a cascade structure of the turbulent velocity field, it is natural to assume that the 
statistical average of velocity differences <5v(r) = v(x + r) — v(x) over a distance |r| follows 
scaling laws 

(|v(x + r) -v(x)| p ) ~r c * (1) 

in the limit of high Reynolds numbers. By taking velocity differences over a distance r, one 
probes objects of corresponding size. 

In addition to this assumption of self-similarity, Kolmogorov [J] also made the seemingly 
intuitive assumption that the local statistics of the velocity field should be independent of 
large-scale flow features, from which it is widely separated in scale. Because the turbulent 
state is maintained by a mean energy flux e, the only local scales available are the length r 
and e itself, which leads to the estimate Sv ~ (er) 1 / 3 or 

(^=p/3 . (2) 

At the same time, one obtains an estimate for the Kolmogorov length 

v = (^ 3 A) 1/4 (3) 

where viscosity is important. However, it was only appreciated later [[J that in turbulence 
long-range correlations always exist in spite of local coupling. Namely, large-scale fluctua- 
tions of the velocity field will result in a fluctuating energy transfer, which drives smaller 
scales. As a result, the statistics of the small-scale velocity fluctuations will be influenced 
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by the energy transfer and fluctuations on widely separated scales are correlated, violating 
the fundamental assumption implicit in (Q) and (0). 

Indeed, Kolmogorov || and Obukhov [Q later proposed the existence of corrections to 
the scaling exponents (|2]), 

C P = p/3 + 5C P ,5C P 7^0 (4) 

which were subsequently confirmed experimentally On one hand, careful laboratory 

experiments have been performed at ever higher Reynolds numbers . On the other hand, 
a new method of data analysis |6],|7j has been successful in eliminating part of the effects 
of viscosity. In particular, for the highest moments up to p = 18 significant corrections 
to classical scaling were found, a currently accepted value for the so-called intermittency 
parameter /i being || 

^=-^ = 0.2 , (5) 

which is a 10 % correction. The existence of corrections like (|5p implies that on small scales 
large fluctuations are much more likely to occur than predicted by classical theory. 

This "intermittent" behavior is thus most noticeable in derivatives of the velocity field 
such as the local rate of energy dissipation 

e(x,t) = -- (dui/dxk + du k /dxif . 

Much of the research in turbulence has been devoted to the study of the spatial structure of 
e(x, t) but which will not be considered here. The statistical average of this quantity 

is what we simply called e before. Owing to energy conservation, it must be equal to the 
mean energy transfer. 

The local coupling structure of turbulence has inspired the study of so-called shell models, 
where each octave in wavenumber is represented by a constant number of modes, which are 
only locally coupled. This allows to focus on the implications of local coupling for intermit- 
tent fluctuations, disregarding effects of convection and mixing. The mode representation of 
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a single shell serves as a simple model for the "coherent structures" a turbulent velocity field 
is composed of, and which to date have only been poorly characterized, both experimentally 
and theoretically. 

We caution that this leaves out two important aspects of turbulence, both of which have 
recently been proposed to lie at the heart of intermittent behavior. First, we have assumed 
that coherent structure are simple in the sense that they only possess a single scale. But 
experimental fill , numerical j0|, and theoretical |13| evidence points to the importance of 



long and skinny threads of vorticity. Although their real significance to turbulence is not 
without dispute ||14|| , they have led to quantitative predictions of intermittency exponents 



15], in excellent agreement with experiment. Second, nonlocal interactions are disregarded 



in shell models. On the other hand it has recently been proposed |16| in the framework of 



perturbation theory that non-local interactions can in fact be re-summed to yield correction 
exponents. 

Given the complexity of the turbulence problem, this makes it all the more interesting to 
carefully assess the possibilities for intermittent fluctuations in the case where both effects 
are eliminated. Thus the aim of this paper is to combine previous numerical [0,0] an d 
analytical |BJ|2IJ efforts to gain insight into intermittent fluctuations in models with local 



coupling alone. We focus on a particular class of shell models, introduced in ||17|| , which 
establishes a direct connection with the Navier-Stokes equation. To this end the Navier- 
Stokes equation is projected on a self-similar selection of Fourier modes, which enforce local 
coupling. We will adopt the name REduced Wave vector set Approximation (REWA) here 

0- 

In the present paper, most of our effort is devoted to separating inertial range fluctuations 
from other effects relating to the fact that scale invariance is broken either by an external 
length scale L or by the Kolmogorov length rj. To this end we simulate extremely long 
cascades, covering up to 5 decades in scale. To eliminate viscous effects as much as possible, 
we use the method of extended self-similarity (ESS) |||7J|. In addition, to make sure our 
results are independent of the method of data analysis, we develop a new subgrid model 



using fluctuating eddy viscosities ||21|| . The results are in excellent agreement with those 
found from ESS. Finally, and perhaps most importantly, we use an analytical calculation 
PDfl to predict the correction exponent without adjustable parameters. Again, we find the 
same values within error bars. We show that the stochastic model introduced in []19| gives an 
excellent description of a REWA cascade, both in an equilibrium and non-equilibrium state. 
We compare the equilibrium properties to adjust a single free parameter in the stochastic 
model. 

Thus we reach two objectives: First, we gain analytical insight into the origin of inter- 
mittent fluctuations in a local cascade. Second, we demonstrate for a simple example how 
equilibrium information about the interactions of Fourier modes can be used to compute 
intermittency exponents. 



Unlike shell models with only one complex mode per shell [23-25] exponent corrections 



5( p for REWA cascades are quite small |17|JT8| . Careful studies of inertial range fluctuations 



have found them to be significant |18[], but their numerical value is only about 1/10 of the 
numbers found experimentally []5| . This has lead to the idea that perhaps the experimentally 
observed exponent corrections are not genuine inertial range properties, but result from 



extrapolations of into regimes where stirring or viscous effects are important [17 JIB, 215 

We will not follow up on this idea here. However, we note that the REWA cascade 
misses an essential feature of three-dimensional turbulence and thus can hardly be expected 
to yield quantitative predictions. Namely, in real turbulence, the number of modes within a 
shell proliferates like r~ 2 as one goes to smaller scales, while in cascade models this number 
is constant. This can be remedied by allowing a particular shell to branch out into eight 



sub-shells, which represent eddies of half the original size. It has been argued [0,0 that 
the competition between eddies of the same size is responsible for the much larger growth 
of fluctuations observed in three-dimensional turbulence. The difficulty with this approach 
is that one also has to take convection of spatially localized structures into account. Also, it 
is not obvious how to disentangle interactions which are local in fc-space from those local in 
real space, as to systematically reduce the coupling of the Navier-Stokes equation to a tree 



of turbulence elements. Recently developed wavelet methods are a promising step , but 
they have so far been used only for data analysis |29|] . 

The paper is organized as follows: In the next section we introduce both the mode- 
reduced approximations of the Navier-Stokes equation and the corresponding Langevin 
models. The inertial range properties of the latter have only one adjustable parameter, 
as explained in the third section. This parameter is determined for a given selection of 
Fourier modes by considering the equilibrium fluctuations of both models. We are thus able 
to predict intermittency exponents of the turbulent state without adjustable parameters. 
The result is compared with numerical simulations of the REWA cascades in the fourth 
section. Exponents are determined by carefully examining various sources of error, and are 
in excellent agreement with the theoretical prediction. In the fifth section, we investigate 
temporal correlations in REWA models to inquire further in the origin of intermittency in 
models with local coupling. This also sheds light on the reasons for the success of the simple 
Langevin model used by us. 
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II. TWO CASCADE MODELS 



A. Reduced wave vector set approximation 

The REWA model ]l7|,[nj is based on the full Fourier-transformed Navier-Stokes equa- 
tion within a volume of periodicity (27rL) 3 . In order to restrict the excited Fourier-modes of 
the turbulent velocity field to a numerically tractable number, the Navier-Stokes equation 
is projected onto a self-similar set of wave vectors /C = ^Q- Each of the wave vector shells 
/Q represents an octave of wave numbers. The shell /Co describes the turbulent motion of 
the large eddies which are of the order of the outer length scale L. This shell is defined 
by N wave vectors Iq . /Co = {kj : i = 1,...,N}. Starting with the generating shell 
/C , the other shells /Q are found by a successive rescaling of /C with a scaling factor 2: 
/Q = 2^/Co- Thus each /Q consists of the iV scaled wave vectors 2 £ k J -°' ) , i = 1, . . . , N. The 
shell JCi represents eddies at length scales r ~ 2~ e L, i.e. to smaller and smaller eddies as 
the shell index I increases. At scales r ~ r\ the fluid motion is damped by viscosity z/, thus 
preventing the generation of infinitely small scales. Hence we only need to simulate shells 
JCe,£ < £ u , where t v ~ log 2 (L/?7) is chosen such that the amplitudes in /C^ are effectively 
zero. In this representation the Navier-Stokes equation for incompressible fluids reads for 
alike /C = Uto^ : 
d 

— Ui(k, t) = -iM ijk (k) u j(Pi ^)Mq ; *) - vk 2 Ui(k, t) + /-(k, t) (6a) 

^ p.qec 

k=p+q 

k-u(k,t) = . (6b) 

The coupling tensor Mij k (k) = [kjPikik) + k k Pij(k)] /2 with the projector P»fc(k) = 5 ik — 
kik k /k 2 is symmetric in j, k and Mij k (k) = — Mij k (— k). The inertial part of equation fl6a]) 
consists of all triadic interactions between modes with k = p + q. They are the same as 
in the full Navier-Stokes equation for this triad. The velocity field is driven by an external 
force f(k, t) which simulates the energy input through a large-scale instability. 
Within this approximation scheme the energy of a shell is 
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W) = \ E l u (MI 2 > 

z ke/Q 



(7) 



and in the absence of any viscous or external driving the total energy of the flow field 
Etot{t) = Yd=QEi(t) is conserved. The choice of generating wave vectors kf^ determines 
the possible triad interactions. This choice must at least guarantee energy transfer between 



shells and some mixing within a shell. In |rT|,I8| different choices for wavenumber sets /Co 
are investigated. The larger the number N of wave numbers, the more effective the energy 
transfer. Usually one selects directions in k-space to be distributed evenly over a sphere. 
However, there are different possibilities which change the relative importance of intra-shell 
versus inter-shell couplings. In this paper, we are going to investigate two different wave 
vector sets, with N = 2Q and N = 74, which we call the small and the large wave vector set, 
respectively. In Fig. 1 a two-dimensional projection of both sets is plotted. The large wave 
vector set also contains some next-to-nearest neighbor interactions between shells, which we 
put to zero here, since they contribute little to the energy transfer. The small set allows 120 
different interacting triads, the large set 501 triads, 333 of which are between shells. 
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FIG. 1. A two-dimensional projection of the /c-vectors in shell /Co for both the REWA models 
considered here. The small set (x) contains all vectors with -1, 0, and 1 as components. The large 
set (0)> i n addition, contains combinations with ±1/2 and ±2. 



Since in the models we consider energy transfer is purely local, the shell energies 
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Ei(t), £ = 0, . . . ,£ u only change in response to energy influx Tg_i_>£ from above and en- 
ergy outfhix T(^( + i to the lower shell. In addition, there is a rate of viscous dissipation 
Tg dlss \t) which is concentrated on small scales, and a rate of energy input T(] n (t), which 
feeds the top level only, cf. Fig.^j. 




FIG. 2. The structure of a local cascade. Eddies of size r ~ 2 L are represented by their total 
energy Eg. Only modes of neighboring shells interact, leading to a local energy transfer T^^ + i{t). 
The cascade is driven by injecting energy into the largest scale with rate Tq \t). The turbulent 
motion is damped by viscous dissipation at a rate T^ d%ss \t). 



From equation (|6a|) we find an energy balance equation which governs the time evolution 
of the shell energies Eg{t) 

jE.it) = T M (t) - T^ e+1 (t) + if ss \t) + Tt\t)S m . (8) 

The different transfer terms are found to be 

(t)=2i M ijk (k)u*(k,t)uj(p,t)u k (q,t) (9a) 

A («+l) 

(0 

Tt\t)= £ Re(u*(k,t)-f(k,t)) (9b) 
ke/Co 

Tf ss \t) = -v k 2 \u(k,t)\ 2 . (9c) 

k€/Q 

In equation (Pal) J] A ('+i) indicates the summation over all next-neighbor triads k = p + q 

(O 

with k e JCg, p e /Q+i and q e Kg U JCe+i- 
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The driving force f(k, t) is assumed to act only on the largest scales, and controls the 
rate of energy input TJj (i). As in Reference [[nj we choose f (k, t) to ensure constant energy 
input Tq 1 ^ = e : 

i{Kt) = 2NHKtW forallke/c ° (10a) 

f(k,t) = for all k £ Kq . (10b) 

This leads to a stationary cascade whose statistical properties are governed by the com- 
plicated chaotic dynamics of the nonlinear mode interactions. Owing to energy conservation, 
viscous dissipation equals the energy input on average. The Reynolds number is given by 

LU eL 2 , . 

since T = (E ) /e sets a typical turnover time scale of the energy on the highest level. We 
believe T to be of particular relevance, since the large-scale fluctuations of the energy will 



turn out to be responsible for the intermittent behavior we are interested in. In [18], for 
example, time is measured in units of L 2//3 e _1//3 , which typically comes out to be l/10th of 
the turnover time of the energy T. As we are going to see below, this is rather a measure 
of the turnover times of the individual Fourier modes. Figure ^| shows the scaling of the 
mean energy in a log-log plot at a Reynolds number of 4.2 • 10 5 . The inertial range extends 
over three decades, where a power law very close to the prediction of classical scaling is 
seen. Below the 10th level the energies drop sharply due to viscous dissipation. In Section 



1 1 J we are going to turn our attention to the small corrections to 2/3-scaling, hardly visible 
in Fig. [3|. Still, there are considerable fluctuations in this model, as evidenced by the plot 
of the energy transfer in Fig. £|. Typical excursions from the average, which is normalized 
to one, are quite large. Ultimately, these fluctuations are responsible for the intermittency 
corrections we are going to observe. 
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FIG. 3. Log-log-plot of the structure function = {En) versus level number for the small 
cascade. At large scales, where the influence of dissipation is negligible, classical scaling is observed. 
At small scales the turbulent motion is damped by viscosity. The Reynolds number is Re = 4.2 -10 5 . 




1.0 2.0 3.0 4.0 

t/T 

FIG. 4. Time evolution of the energy transfer To_>i(£) from shell /Co to shell K-\. Typical 
excursions are large compared with the mean value (To_>i) = e. The time is given in units of 
T = (Eq) /e, where (Eq) is the mean energy of the top level. 



Thus within the REWA-cascade we are able to numerically analyze the influence of 
fluctuations on the stationary statistical properties of a cascade with local energy transfer 
on the basis of the Navier-Stokes equation. In Fig. § we plot the time evolution of the 
energy on the second level of the cascade. One observes 
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FIG. 5. Energy of a shell (£ = 2) as a function of time. The rapid fluctuations come from 
the motion of individual Fourier modes. A much longer time scale is revealed by performing a 
floating average over one turnover time of the second level (bold line). The time is given in units 
of T = {E ) /e. 



short-scale fluctuations, which result from the motion of individual Fourier modes within 
one cascade level. However, performing a floating average reveals a second time scale, which 
is of the same order as the turnover time of the top level. As we are going to see in Section 



FVl , this disparity of time scales is even more pronounced on lower levels. The physics idea 
is the same as in the microscopic foundation of hydrodynamics, where conserved quantities 
are assumed to move on much slower time scales than individual particles. This motivates 
us to consider the energy as the only dynamical variable of each shell, and to represent the 
rapid fluctuations of Fig. |5| by a white-noise Langevin force. In this approximation we still 
hope to capture the rare, large-scale events characteristic of intermittent fluctuations, since 
the conserved quantity is the "slow" variable of the system. Similar ideas have also been 
advanced for the conservative dynamics of a non-equilibrium statistical mechanical system 
|3TJ. 
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B. The Langevin cascade 



In this model we take a phenomenological view of the process of energy transfer. The 
chaotic dynamics of the REWA cascade is modeled by a stochastic equation. We make 
sure to include the main physical features of energy conservation and local coupling. In 
particular, the dynamics is simple enough to allow for analytical insight into the effects of 
fluctuating energy transfer |20| . 

As in the preceding REWA-cascade, the turbulent flow field is described by a sequence of 
eddies decaying successively (Figure 0). The eddies at length scales r ~ 2~ e L are represented 
by their energy Ee(t). As before we restrict ourselves to local energy transfer, and thus the 
time evolution of the shell energies E^{t) is governed by equation The crucial step is 
of course to choose an appropriate energy transfer 7}_^+i(t). For simplicity, we restrict 
ourselves to a Langevin process with a white noise force. Thus the local transfer T^ +1 (i) 
is split into a deterministic and a stochastic part T^ + i(t) = T^f +1 (t) + T^° e ^(t) where 
both parts should depend only on the local length scale 2~ e L and the neighboring energies 
Eg and Et+\. The most general form dimensionally consistent with this has been given in 
20| . For simplicity, here we restrict ourselves to the specific form 



TiH.it) = D^j- (Ef(t) - E%\(t)) (12a) 

TtXKt) = R E 5 / 8 (t)Et +1 (tf/% +1 (t) (12b) 

Tt\t)=e5 m (12c) 
Tf is \t) = -u{2-"L)- 2 E, . (12d) 

The white noise is represented by i.e. (£t(t)) = and (£t{t) = 25w5(t — t'). We 
use Ito's |H] definition in equation ( |12b| ). To understand the dimensions appearing in (|12|). 



note that ui ~ E 1 ^ 2 is a local velocity scale and k ~ 2 e /L is a wavenumber. Thus ( |1 2a| ) 



dimensionally represents the energy transfer (|9a|) . In ( |12b| ) the powers are different, since 
£ carries an additional dimension of l/time 1//2 . It follows from ( |12a| ) that the sign of the 
deterministic energy transfer depends on which of the neighboring energies Eg_ or E^ + i are 
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greater. If for example Eg is larger, T^^ +1 (t) is positive, depleting Eg in favor of Eg + i. Hence 
the deterministic part tends to equilibrate the energy among the shells. The stochastic part, 
on the other hand, is symmetric with respect to the two levels i and i + 1. This reflects our 
expectation that in equilibrium it is equally probable for energy to be scattered up or down 
the cascade. 

The combined effect of ( |12a| ) and ( |12b| ) is that without driving, energies fluctuate around 
a common mean value. This equipartition of energy in equilibrium is precisely what has 
been predicted on the basis of the Navier-Stokes equation |32l , [33| . The only free parameters 
appearing in the transfer are thus the amplitudes D and R. If R is put to zero, the motion 
is purely deterministic, and one obtains the simple solution 

= C2-<W with C = (^) 2/3 • (13) 

This corresponds to a classical Kolmogorov solution with no fluctuations in the transfer. The 
amplitude D of the deterministic part is a measure of the effectiveness of energy transfer. 



On the other hand R measures the size of fluctuations. In p(J it is shown that a finite R 
necessarily leads to intermittency corrections in the exponents. In the next section we are 
going to determine the model parameters for the two REWA cascades we are considering. 
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III. DETERMINATION OF MODEL PARAMETERS 



The aim of this section is to explain the significance of the parameters appearing in the 
Langevin model. We show that only the combination R/D 1 ^ 2 determines the nonequilibrium 
fluctuations in the inertial range. But it is the same combination which also sets the width of 
the equilibrium distribution of energies, if the chain of shells is not driven. Thus we are able 
to fix all the parameters necessary to describe the nonequilibrium state solely by measuring 
equilibrium properties. 

The physical parameters of the turbulent cascade are the length scale L of the highest 
level, the rate of energy input e, and the viscosity v. The properties of the energy transfer are 
determined by the dimensionless strength of the deterministic part D and of the stochastic 
part R. In the previous section we have seen that = (2eL/ D) 2 ^ 3 sets an energy scale 
for the highest level, and thus T = /e is a time scale. Both scales can be used to 
non-dimensionalize (|T2]) , giving 

f^Ut) = 2 m (E'fit) - (14a) 
TtXht) = -^2^/ 2 ^ /8 (t)^ + i(t) 5/8 6 + i(t) (14b) 
f (m) (t) = l (14c) 
f l {diss \t) = -2 2 \Re)' 1 E l . (14d) 

(14e) 

Hence R/D lj/2 is the only parameter characterizing the inertial range dynamics, which can 
be understood as follows: At any level £, a time scale of deterministic transport is set by 
td ~ Ei/e ~ E e 1 ^ 2 r/D, where Ei is a typical energy scale and r = 2~ e L. This deterministic 
transport competes with diffusion of energy introduced by the stochastic part of the energy 
transfer. Namely a time scale over which the energy Ee can diffuse away is given by tr ~ 
(Ee/A) 2 , where A = Rer^^E^ 4 is the amplitude of the noise term. Plugging in the 
expression from ( |12b|) one ends up with ~ rEJ 1 ^ 2 /R 2 , and hence 

tr D_ 

T D ~ R 2 ' 
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In particular, the relative importance of deterministic transport and diffusion is independent 
of the level, and only depends on the combination of D and R given above. 

We are now left to determine Rj D 1 ! 2 from an experiment which is independent of the 
nonequilibrium state. To that end we consider a long chain of shell elements, for which 
energy input as well as viscous dissipation has been turned off (e = v — 0). As a result, the 
energy will perform fluctuations around some mean value (Eg). Using a similar argument 
as above, the probability distribution p(Eg/ (Ei)) for the shell energies will only depend on 
R/D 1 / 2 for the Langevin model. Figure |6| shows the probability distribution for one level 
of a cascade with 7 shells. Except for some end effects at the lowest level, the distribution 
turns out to be level- independent. The large REWA cascade is compared 

10* 
10 1 
10° 

itf 1 

10" 2 
10 3 

0.5 0.75 1.0 1.25 1.5 

E /<E ) 

FIG. 6. Equilibrium distribution of the energy in the top level of the large REWA-cascade 
(N = 74), consisting of 7 shells. By adjusting the only parameter R/D 1 ^ 2 , the distribution of the 
Langevin cascade can be made to agree within statistical error. The results for R/D 1 / 2 for both 
wave vector sets are given in (|lq) . 

with a simulation of the Langevin cascade. Once R/D 1 / 2 is adjusted, there is an almost 
perfect match between the two models. This shows that the simple stochastic dynamics 
proposed here models the equilibrium distribution of the chaotic fluctuations of Fourier 
modes very well. Most importantly, we have determined the only adjustable parameter of 
the energy transfer of the Langevin cascade corresponding to the two model systems we are 
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considering: 

R 

— — = 1.729 ■ KT 1 for the small cascade (N = 26), (15a) 
R 

— — r = 9.968 • 1(T 2 for the large cascade (N = 74) . (15b) 
Jj 1 / 2 

As seen from (plj), the REWA cascade with a larger number of modes has smaller fluctua- 
tions. This is not surprising, since large excursions of the energy correspond to a coherent 
motion of the individual Fourier modes. The larger the number of modes, the harder this is 
to achieve, since the random motion of individual modes tends to destroy correlations. 
In the next section, the fit ( |15a| ),( |IT)T5|) will be used to compare nonequilibrium properties 



in the inertial range. All information about possible differences between the small and the 
large cascade has been condensed into a single number. 
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IV. INTERMITTENCY CORRECTIONS 



We now turn to the inertial range fluctuations of the two cascade models. To allow for a 
direct comparison, we focus on the scaling of the energies Eg of one shell. The analogue of 
the moments of the velocity field usually considered in turbulence are the structure functions 
based on the energy 

= (E p t /2 ) ~ 2-^ 1 . (16) 

We concentrate on the small corrections S( p to the scaling exponents as predicted on di- 
mensional grounds. As seen in Fig. these corrections are extremely small, so that no 
significant deviation from classical scaling is seen on the scale of the figure. 



In order to know what to expect, we use the result of Ref. [20] where we have computed 
the scaling behavior of a general class of stochastic models in a perturbation expansion. To 
lowest order, the exponent corrections are given by the quadratic dependence 

*Cp = -^(p-3) • (17) 



Plugging the specific form of the energy transfer (|14a ), ( |14b| ) into the formulae given in |20 
we find 



R 



2 



/i=0 - 4 H^J • (18) 

For the small cascade this leads to /x = 0.013, which is about 1/10 of the experimental value 
accepted for three-dimensional turbulence || . Therefore, it is more than ordinarily difficult 
to measure the scaling exponents with sufficient accuracy to obtain significant answers for 
the deviations from classical scaling. On the other hand, the number of modes in a REWA 
cascade being greatly reduced as compared with the full Navier-Stokes equation, we are able 
to simulate cascades with up to 17 levels, corresponding to 5 decades in scale. 

To obtain reliable results within the accuracy needed, it is essential to disentangle statis- 
tical errors from systematic errors, introduced through finite-size effects or viscous damping. 
To that end we individually assign statistical errors to every average taken. Figure |7| 
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FIG. 7. Convergence of the structure function = on the 11th step of the small 

REWA cascade. We show the evolution of D\ as function of the averaging time T av in comparison 
with the stationary values £)( 18 ' stat ^ ( w hich are calculated by integrating over all initial values in the 
data set. The error bars give the statistical error of the data as estimated by (|T9|). The averaging 
time T av is given in units of T = (Eq) /e. 




shows the convergence of the 9th moment of the energy on the 11th step of the small 
REWA cascade, which is the lowest level relevant to our fits. We plot the mean value, 
averaged over the time given. Evidently the fluctuations of this mean value get smaller as 
the averaging time T av is increased. To obtain a quantitative measure of the uncertainty 
of the T-average of some quantity x, we also consider the ensemble of averages obtained 
with different initial conditions. The result for the mean value and variance of this ensemble 
average is plotted as circles with error bars in Fig. [7], where E\ takes the place of x. The 



variance is found equivalently as an integral over the temporal correlation function [34 



1 rT av rT av r 2 1 

a x (T av ) = — / / (x(ti)x(t 2 )) - {x) dTidr 2 . (19) 

T av JO JO 1 J 

This variance cr(T av ) is seen to give a reasonable approximation to the fluctuations of the 
temporal average. The variance for the largest averaging time available has been taken as 
the statistical error of a measured temporal average. 

Next we consider systematic errors in the computation of the scaling exponents. Devi- 
ations from power law scaling are expected to occur on both ends of the cascade and have 
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been studied extensively []17|p5|,p!8[ . First, fluctuations on the highest level are suppressed, 
since there is no coupling to a higher level, but deterministic energy input instead. Second, 
in the viscous subrange the energy is increasingly depleted by viscosity, leading to even more 
drastic effects on the spectra. Both effects are strongest for the highest moments, which are 
most sensitive to large fluctuations. The 18th order structure function for a typical run of 
the small cascade is plotted as diamonds in Fig. || The average on the highest level is 
considerably lower than expected from scaling, as fluctuations are suppressed. 
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(O) versus level number (Standard) in comparison with extended self similarity (ESS) scaling 

~ (18) ~ (3) ~ (^) 

of Dy°> (□) versus D), /Dq . The structure functions are normalized by their power law fits: 
D < f> norm) = D^/(b^2-^) in the standard plot and D^' norm) = /{A^ {bf /b®)M&) in 
the ESS plot. Refer to Table | for the values of the exponents. In the ESS plot the scaling range 
is more than doubled. The numerical calculation is performed for the small REWA cascade. 



The cascade only gradually recovers from this suppression, which leads to a faster rise 
in the level of fluctuations and thus to a decrease in the local scaling exponent as reported 
18 1 . Note that we show a "scatter-plot" of the data around the power law D^' ~ 2~^ p , 



111 



which represents our best fit. Thus deviations from this power law are hugely exaggerated 
and would not be visible on a customary logdog-plot. The absolute variation of D\ over 
the range of the plot is 21 orders of magnitude. Owing to the extreme sensitivity of our plot, 
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TABLE I. Correction exponents of the structure functions DY' for p = 2, 6, 12 and 18 for the 
small REWA cascade. The exponents in column ESS are determined by plotting versus 
as suggested by extended self-similarity. To obtain the values in column EV, viscous effects have 
been removed entirely by introducing a fluctuating eddy viscosity. In the last column the prediction 
based on the Langevin model is given. 

viscous damping is visible below the fourth level, even for a cascade with 17 shells, as seen 
in Fig. [| This restricts the inertial range to four levels. We employ two different strategies 
to improve on this situation: (i) We plot the higher order structure functions against the 
third order structure function, as suggested by the ESS method, (ii) We eliminate viscosity 
altogether by introducing a fluctuating eddy viscosity instead. 

It is known P,|7J that extended self-similarity leads to a very considerable improvement 
of the scaling of experimental data. Without any viscous corrections, the scaling behavior 
is unaffected, as seen in Fig. || in the fitting range marked "standard". Below the forth 
level, however, the third order structure function suffers viscous damping quite similar to 
that affecting the other structure functions. Thus by plotting D\ on the abscissa, both 
effects are hoped to largely cancel each other. This is indeed true for the REWA cascade as 
well, where the scaling range has been extended to almost three decades. In Table | we have 
compiled various scaling exponents for the small cascade. Except for the smallest moment, 
highly significant corrections to the scaling exponents are measured. The errors are based 



on a least squares fit with weighted averages [36], based on the statistical errors as explained 
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above. To make sure the results do not depend on our choice of the inertial range, we also 
made fits in ranges other than the one marked "ESS" in Fig. [8[ Namely, we variously 
shifted the fitting range by one level up or down, or took an additional level into account 
on either end. In each case, the values of the exponents were within the errors given in 
Table |. Typical averaging times were 100 turnover times of the largest scale. This is about 
10 times as long as in fT8|, as we base the turnover time on the time scale of the energy. 
For comparison, we also supply the value of the exponent correction as calculated from (|17D 
and on the basis of the known value of R/D 1 ^ 2 , given in (|15|) . There is no adjustable 
parameter in this comparison with the theoretical prediction, since the noise strength was 
determined solely from the equilibrium fluctuations of the REWA cascade. The excellent 
agreement found for all exponents makes us confident that finite size corrections and viscous 
effects have been successfully eliminated, and we are measuring genuine properties of the 
inertial range. 

However, since there is no theory demonstrating that the ESS method leads to the 
correct inertial range scaling behavior, we have nevertheless attempted to eliminate viscous 
effects by a second, and completely independent method. This was done by putting v = 0, 
and instead draining energy from the lowest shell JCe max using an eddy viscosity [^,|T^]. In 
order to mimic inertial energy transfer into the subgrid shells K.£,£ > £ m ax, we preserve the 
coupling structure of the equations as much as possible. We add a term 



D(k,t) 



-d(k,t)2 £ — |u(k,t)|u(k,t) : k6/C w 

(20) 

: k g JC {mit) 



to the inviscid Navier-Stokes equation (|6a|) . The cut-off shell K,( ut) C fcg max contains all the 
wave vectors of Ki max which interact directly with modes of the shell which is not resolved. 
The problem of this procedure is that (i) the fluctuations in the subgrid scales are not 
accounted for, and (ii) the amplitude to choose for d(k, t) is not known. Both problems are 
addressed by using a method inspired by the work in |21|. The idea is to adjust d(k, t) such 
that 
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u(k)| : 



^(lu(k')l 3 



for k e /Q„ 



k' = -k 



(21) 



which is valid exactly in a perfectly scale- invariant cascade [|T7| , analogous to the Kolmogorov 
structure equation [|1|]. The wave vector k' e K<t max -i in (0) is a scaled copy of k e JCe max . 
In the long time limit, each d(k, t) will converge to some average value, which is determined 
self-consistently from the dynamics of the cascade. But in that case the fluctuations of the 
subgrid scales would be missing. Therefore, we instead took the averages in fl2~ID over just 10 
turnover times of the last resolved shell JCe max . When (|u(k)| 3 ) / (|u(k')| 3 ) was larger than 
1/2, we increased d(k,t) by 5 %, otherwise d(k,t) was decreased by the same percentage. 
Thus all the d(k, t) reflect the fluctuations occurring at the end of the cascade, as seen in Fig. 
|9] for one of the amplitudes of the small cascade. The average of d(k, t) becomes stationary 
in the long time limit, but fluctuations are considerable, as expected in an intermittent 
cascade. 




FIG. 9. Time evolution of the amplitude d(k, t) of the eddy damping (]2^) for k = 2 lmax (l, 1, 1). 
The fluctuations mirror the highly intermittent behavior of the velocity modes on the 14th level. 
The numerical calculation is performed with the small wave vector set. Time is given in units of 
T = (Eq) jt. 
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FIG. 10. Anomalous scaling of the structure function D\ , p = 18 for the small REWA 
cascade, with the eddy viscosity (|20|) on the lowest level. We show the scatter of the data around 
the fitted power law by plotting f)^ norm ^ = ^ 18) /(& (18) 2 _Cl8 ^)- Classical Kolmo gorov scaling is 
shown for comparison. 



Figure [10] shows the 9th moments of the energy as in Fig. |8|, but for the eddy-damped 
cascade. End effects at the small scales are very small, less substantial than on the largest 
scales. Also shown are the error bars resulting from the statistical estimate described before. 
The result of the fit over the scaling range indicated is given in Table | for various moments. 
Again, the range of our fit was also varied, and the results are consistent with the errors 
given. The values of the exponent corrections are in excellent agreement with the values 
obtained from ESS. Furthermore, the error is even less than before, giving very significant 
deviations from classical scaling. This is also seen from the dashed line in Fig. [TO, which 
represents classical Kolmogorov scaling. Finally it should be noted that even for the long 
averaging times we use, some imbalances in the cascade remain, which make D\ deviate 
from its exact power law behavior of 2~ e . These deviations decrease in time and turn out not 
to be significant for many of our runs. However, we decided to normalize our results so as 
to substract remaining imbalances. So strictly speaking the values given for the correction 
exponents in Table | are 

Cp Cs 3 ' 
24 



REWA: 


N = 


= 26 




1.3 


io- 


" 2 ± 1 


10~ 3 


Langevin: 


R/D 1 / 2 = 


= 1.729 


lCT 1 


1.3 


10~ 


-2 




REWA: 


N = 


= 74 




5.3 


10- 


" 3 ± 1 


10~ 3 


Langevin: 


R/D 1 / 2 = 


= 9.968 


1(T 2 


4.2 


10' 


-3 





TABLE II. Comparison of our best numerical estimate of of the intermittency correction 
fi = —5^6 for the REWA cascades with the theoretical prediction. In both cases, with the predicted 
exponents differing by a factor of three, agreement is within error bars. 



The same procedure has been adopted in previous work on the subject |I8| 



Since the best error estimates are obtained by using an eddy-damped cascade, we are 
using this method for our final comparison between numerical simulations and the theoretical 
prediction of the Langevin model. The result of this comparison for /x = —5(q is found in 
Table [Tl] for the small and the large cascade. In both cases, numerics of the REWA cascade 
and theory agrees within error bars, the size of the exponent corrections differing by more 
than a factor of three between the small and the large cascade. This underscores 
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FIG. 11. Intermittency corrections for the two REWA cascades in comparison with the predic- 
tion of theory. The error bars of the numerical values are based on a weighted least square fit, and 
are given exemplary for the highest moment. 
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the consistency of our results and demonstrates that the Langevin model captures all 
the essential physics responsible for the build-up of fluctuations in a local cascade. The 
same message is contained in Fig. [TI|, which summarizes the exponent corrections for both 
models. A better understanding of why a model, which has the energy as its only mode 
works so well, is supplied by a study of the temporal correlations. At the same time it gives 
insight into the origin of intermittency itself, at least in models with local coupling. 
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V. TEMPORAL CORRELATIONS 



In the previous sections we have looked only at equal time correlations. Although they 
contain information about the fluctuations of the cascade, their information about the dy- 
namics responsible for these fluctuations is very indirect. So ultimately one has to look at 
temporal correlations as well to understand the dynamical origins of intermittency. In the 
case of the REWA cascade, one can look at the fluctuations of the individual Fourier modes, 

C^(k,t) = (u*(k,0)- U (k,t)), (22) 

just as one does in studies of the full Navier-Stokes equation. We also expect the temporal 
correlations of the total energy within a cascade step to be of particular significance, 

Ci B) (t) = (E l (0)E t (t))-(E t ) 3 , (23) 

since the conservation properties of the energy are responsible for maintaining a turbulent 
state. 

The standard guess for the scale dependence of correlation times is based on the Kol- 
mogorov picture. Namely, assuming that the typical correlation time is a local quantity, the 
only combination of the length scale 2~ e L and the mean energy transfer e having dimensions 
of time leads to 



n 



2- l Lf/t 1/3 . (24) 



We test this idea by computing the temporal correlation ( f22f ) of a particular mode of the 
REWA cascade, for three different shells, as plotted in Fig. O. 
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FIG. 12. The real part of the autocorrelation function C^(k,t) = (u*(k, 0) • u(k, t)) normal- 
ized to its equal time value CW(k, 0) for the mode k = 2^(1, 1, 1) of the small cascade. Shown are 
the three top levels. The correlation function decays very quickly, reflecting the chaotic behavior 
of the velocity field. The time, given in units of T = (Eq) /e, is rescaled according to the classical 
prediction ~ 2 -2 ^/ 3 for the decay time. 



The abscissa is rescaled according to fl24|) , which makes all three correlation functions 



fall atop of each other. So it seems that if corrections to (p4 ) exist, they are so small that 
they need more sophisticated methods to be revealed. 

Next, we look at the temporal correlations of the energy on the highest level, plotted in 
Fig. [T3|. 
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FIG. 13. The autocorrelation function Cg (t) of the energy on the top level, normalized to its 
equal time value Cq (0) for the small cascade. The energy de-correlates by more than a factor of 
10 more slowly than a single velocity mode. The time is given in units of T = (Eq) je. 




Apart from differences in shape, the remarkable feature is that the typical decay time 
is more than 10 times as long as that of individual Fourier modes. The temporal corre- 
lation is apparently dominated by the long-time fluctuations of the energy already noted 
in Fig. [|. The rapid fluctuations representative of the individual Fourier modes appear 
to be completely uncorrelated. The physical reason is that many different modes within a 
shell contribute to these rapid fluctuations, which are de-correlated through many random 
"collisions" . Just like in a gas of particles, the complicated interaction between many modes 
tends to randomize the individual motions. 

If the small-scale fluctuations were a true random walk, the energy would slowly diffuse 
to take arbitrary values. But eventually the tendency of the dynamics to restore equilibrium 
will drive the energy back to its mean value. This requires a coherent motion of many 
individual Fourier modes, whose correlations need some time to build up. Thus the long 
time scale in the motion of the energy. So far the same argument would apply for a cascade 
in equilibrium, in the absence of driving. This is shown in Fig. [Tl for a cascade of 8 shells, 
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FIG. 14. Correlations times t« u,E ' for the small cascade in equilibrium, as a function of the 
level number. The dashed line corresponds to the theoretical prediction of (|25|). The time scale of 
the energy is longer than that of an individual mode by a factor of 4 on all levels. 



which all perform fluctuations whose static distribution was shown in Fig. []. Since the 
energies of all shells fluctuate around some common value E av , and the local length scale is 
2~ e L, for dimensional reasons the time scale must be 



(e qu) ^ E .y 2 2 _, L (25) 



The correlation times shown in Fig. [14] for both the energy and one of the velocity modes 
is calculated according to 

(E,u) _ Io°Ce E ' U \t) dt 

and are found to corroborate the scaling law (|25|) . As in the case of the top level of a 
Kolmogorov cascade, the correlation time of the energy is slower compared to an individual 
Fourier mode. 

This changes fundamentally when looking at lower shells of a non-equilibrium cascade. 
For a turbulent cascade, one expects to recover the scaling (P4*D, since the correlation of the 
energies are denned in terms of local quantities only. But remarkably, the scaling law ( P4"| ) 
does not even approximately describe the scaling of r £ , which rather follows the power law 

rf ] ~ 2~ a£ , a = 0.094, (27) 
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as seen in Fig. 
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FIG. 15. Scaling of the correlation times t) of the energy correlation function C\ (t) for the 
small cascade. The fit corresponds to Tg ~ 2 -0 ' 094 ^, which is a very small decrease compared to 
the classical prediction of Tg ~ 2 -2 ^/ 3 . 



This means the time scale of the motion of the energy hardly gets shorter at small scales. 
The reason lies within the non-equilibrium properties of the cascade: Once a large fluctuation 
of the energy has built up on the highest level, it can only relax by being transported to 
a lower level. Thus the same long time dependence is imprinted on the lower level, which 
again can only be transferred to the next level. Thus on a given level, all time scales of 
the levels lying above appear, and apart form the single rescaling (]2"? ) the shape of the 
correlation function changes as well. This corresponds directly to the origin of intermittent 
fluctuations itself: fluctuations amplify because lower shells are driven by slowly varying 
energy input of the higher shells. Consequently fluctuations "ride" atop of the long-scale 
fluctuations and amplify. By looking at correlations between individual Fourier modes like 
none of these long-range correlations are revealed. Correlations between the coherent 
structures of a turbulent velocity field will always be masked by the incoherent fluctuations 
of the individual modes. 

It is interesting to note that in a numerical simulation of isotropic turbulence Yeung and 



Pope [22] have found very strong deviation from classical scaling as well. They looked at the 
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Reynolds number dependence of the acceleration variance, which is a Lagrangian quantity. 
Like in our case, including fluctuations of the energy transfer according to Kolomogorov's 
refined similarity hypothesis || cannot account for the corrections found. However, a mean- 
ingful definition of a Lagrangian quantity within our local approximation is difficult, since 
convection cannot be described properly. Therefore we do not know how to relate our 
findings with the results of Ref. directly. 

The second consequence of the extremely slow de-correlation of the energy is that the 
ratio of time scales of the energy and that of individual modes rapidly becomes larger on 
smaller scales. In the Langevin description the white noise term represents the irregular 
motion of individual Fourier modes, while our main interest lies with the large-scale fluctu- 
ations of the energy. In the scaling limit we are interested in the disparity between these 
scales becomes infinitely large. Thus a white-noise description should become better and 
better in the relevant limit. 
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VI. DISCUSSION 



We have investigated mode-reduced approximations of the Navier-Stokes equation and 
found them to have anomalous scaling exponents in the inertial range. Intermittent fluctua- 
tions come about through rare excursions of the energy from its mean value, which originate 
from the top levels of the cascade. They are thus well described by a white-noise Langevin 
process, whose random forcing represents the motion of individual Fourier modes. The an- 
alytical solution shows |21J that such a Langevin cascade necessarily exhibits corrections to 
the classical scaling exponents. It is thus hard to see how any cascade with local coupling 
could avoid intermittency corrections in the inertial range, since for a chaotic motion there 



will always be fluctuations, and thus D/R ^ 0. In fact, as pointed out by Kraichnan |38 



the only mechanism by which such fluctuations could be avoided is by sufficiently strong 
mixing in space, which in /c-space would be represented by non-local interactions. The fact 
that our simple Langevin model gives an excellent description of a complicated network of 
Fourier modes is also highlighted by the excellent agreement in an equilibrium state. We 
have recently extended the comparison between the Langevin and the REWA model to the 
formation of singularities in the absence of viscosity |39 |. In this case, which represents a 
state even farther from equilibrium than a Kolmogorov cascade, the two models still com- 
pare quantitatively. This allows for an analytical description of the singularities of the Euler 
equation in the approximation of local coupling. 

The second central point of our paper is to demonstrate that intermittency exponents 
can be calculated for a complex cascade with nonlinear, chaotic dynamics by analytical 
means. The only piece of information one needs about the turbulent flow concerns the local 
energy transfer. With this information, the global non-equilibrium properties determine the 
intermittency exponents, which is taken care of by the analytical calculation. In a very 



interesting paper, Olla [^| has recently carried out the same program directly from the 
Navier-Stokes equation. He assumes that the broadening of the distribution of the velocity 
field from a scale r to r/2 is adequately described in a Gaussian approximation. Thus he 
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is able to determine the analogue of the noise strength R/ D 1 / 2 from perturbation theory. 
However, a major problem of Olla's work is that he also does not deal adequately with the 
spatial structure of turbulence, but rather maps the problem on the same linear cascade 
structure we are considering here. 

The most pressing problem remaining with the present mode representations is that the 
intermittency effects are quite small. Thus some important physical features are still missing. 
The most likely reason is the fact that a given eddy only feeds a single substructure, instead 
of branching out to several smaller structures overlapping with it in physical space. It has 
been shown |2T| , [T9|j that such a tree structure may strongly enhance intermittent fluctuations, 
owing to competition between different eddies on the same level. The procedures of the 
present paper can be directly carried over to this physical situation, as analytical solutions 



of the Langevin model with tree structure are available as well ||40|| . Unfortunately, there 
is no technique available which would link the Navier- Stokes equation to a cascade coupled 
locally in real and in /c-space in a rational way. Since the value of intermittency exponents 



depends significantly on the properties of the local couplings |3q] , no quantitative prediction 
for real three-dimensional turbulence is possible using the more general tree structure. 
An interesting fact revealed by our analysis is that the value of the correction exponents 



depends considerably on the mode selection. There is a tendency, observed earlier [18 



for the size of fluctuations to decrease with the number of modes per shell. In [|T8|, mode 
systems with up to 86 modes per shell have been considered, which also allow for more 
distant interactions between shells. Although it is interesting to look at the significance 
of more distant interactions, it does not remedy the fundamental problem of the present 
approximation, which is that the number of modes does not proliferate on small scales. In 
terms of interactions, there is still no competition between modes of the same scale. Thus 
no conclusions about the true intermittency exponents can be drawn from including more 
modes, as long as the spatial structure is not respected. In fact for a tree structure the 
number of localized boxes would increase like 2 3£ , each with a constant number of modes 
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N. Comparing that with the total number of Fourier modes (47r/3)(2 3 — 1)2 3£ within a 
shell for three-dimensional turbulence, one obtains the estimate iV « 28, which is close to 
the number = 26 of our small cascade JEJ. Thus the small cascade already represents a 
physically reasonable choice for the number of locally interacting Fourier modes. 

Finally we mention the so-called GOY model which contains only one complex 

mode per shell, and which has attracted considerable attention recently. Although it has no 
spatial structure, the intermittency exponents are comparable to those of real turbulence, 
which seems to contradict our above observations. We believe that the mechanisms at work 
are somewhat different here, since there exist pulses, which run down the cascade. These 
global coherent motions are much more effective in transporting fluctuations to smaller 
scales than in REWA cascades, where such structures are destroyed by intra-shell mixing. 
A separation of time scales between the total energy and the individual modes of a shell 
evidently does not exist. Another fact which points to global structures is the dependence 
of intermittency exponents on the form of viscous dissipation ||45|| . Still, this simple model 
remains interesting in particular in view of the analytical theory which has been developed 
for its inertial range properties P6 |, and which has a number of similarities with the one 
developed for the Langevin model AISLED • 



To conclude, we have investigated in detail how local coupling leads to intermittent fluc- 
tuations. In spite of many other interesting proposals, Landau's and Kolmogorov's original 
arguments remain among the most powerful for their understanding. 



35 



ACKNOWLEDGMENTS 

This work is supported by the Sonderforschungsbereich 237 (Unordnung und grosse Fluk- 
tuationen). 



36 



REFERENCES 



[1] A. N. Kolmogorov, "The local structure of turbulence in incompressible vis- 
cous fluid for very large Reynolds' numbers", C. R. Akad. Nauk SSSR 30, 301 
(1941); "On degeneration of isotropic turbulence in an incompressible viscous fluid", 
C. R. Akad. Nauk SSSR 31, 538 (1941); "Dissipation of energy in the locally isotropic 
turbulence", C. R. Akad. Nauk SSSR 32, 16 (1941). 

[2] L. D. Landau and E. M. Lifshitz Fluid Mechanics, (Pergamon, Oxford, 1959; third 
edition, 1984). 

[3] A. N. Kolmogorov, "A refinement of previous hypotheses concerning the local structure 
of turbulence in a viscous incompressible fluid at high Reynolds number", J. Fluid Mech. 
13, 83 (1962). 

[4] A. M. Obukhov, "Some specific features of atmospheric turbulence" , J. Fluid Mech. 13, 
77 (1962). 

[5] F. Anselmet, Y. Gagne, E. J. Hopfinger and R. Antonia, "High-order structure functions 
in turbulent shear flows", J. Fluid Mech. 140, 63 (1984). 

[6] R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, F. Massaioli and S. Succi, "Extended 
self-similarity in turbulent flows", Phys. Rev. E 48, R29 (1993). 

[7] R. Benzi, S. Ciliberto, C. Baudet, G. R. Chavarria and R. Tripiccione, "Extended self- 
similarity in the dissipation range of fully developed turbulence", Europhys. Lett. 24, 
275 (1993). 

[8] J. Herweijer and W. van de Water, "Universal shape of scaling functions in turbulence", 
Phys. Rev. Lett. 74, 4651 (1995). 

[9] B. Castaing, Y. Gagne and E. J. Hopfinger, "Velocity probability density functions of 
high Reynolds number turbulence", Physica D 46, 177 (1990). 



37 



[10] M. Nelkin, "What do we know about self-similarity in fluid turbulence?" , J. Stat. Phys. 
54, 1 (1989). 

[11] S. Douady, Y. Couder, and M. E. Brachet, "Direct observation of the intermittency of 
intense vorticity filaments in turbulence", Phys. Rev. Lett. 67, 983 (1991). 

[12] S. Kida and K. Ohkitany, "Spatio-temporal intermittency and instability of a forced 
turbulence", Phys. Fluids A 4, 1018 (1992). 

[13] H. K. Moffatt, S. Kida, and K. Ohkitany, "Stretched vortices - the sinews of turbulence; 
large-Reynolds-number asymptotics" , J. Fluid Mech. 259, 241 (1994). 

[14] J. Jimenez, A. A. Wray, P. G. Saffman, and R. S. Rogallo, " The structure of intense 
vorticity in inhomogeneous isotropic turbulence", J. Fluid Mech. 255, 65 (1993). 

[15] Z.-S. She and E. Leveque, "Universal scaling laws in fully developed turbulence", Phys. 
Rev. Lett. 72, 336 (1994). 

[16] V. L'vov and I. Procaccia, "Exact resummations in the theory of hydrodynamic turbu- 
lence. I. The ball of locality and normal scaling", Phys. Rev. E 52, 3840 (1995); "II. A 
ladder to anomalous scaling", Phys. Rev. E 52, 3858 (1995). 

[17] J. Eggers and S. Grossmann, "Does deterministic chaos imply intermittency in fully 
developed turbulence?", Phys. Fluids A 3, 1958 (1991). 

[18] S. Grossmann and D. Lohse, "Scale resolved intermittency in turbulence", Phys. Fluids 
6, 611 (1994). 

[19] J. Eggers, "Intermittency in dynamical models of turbulence", Phys. Rev. A 46, 1951 
(1992). 

[20] J. Eggers, "Multifractal scaling from nonlinear turbulence dynamics: analytical meth- 
ods", Phys. Rev. E 50, 285 (1994). 

[21] J. A. Domaradzki and W. Liu, "Approximation of subgrid-scale energy transfer based 



38 



on the dynamics of resolved scales of turbulence", Phys. Fluids 7, 2025 (1995). 

[22] P. K. Yeung and S. B. Pope, "Lagrangian statistics from direct numerical simulations 
of isotropic turbulence", J. Fluid Mech. 207, 531 (1989). 

[23] A. M. Obukhov, "Some general properties of equations describing the dynamics of the 
atmosphere", Atmos. Ocean. Phys. 7, 471 (1971). 

[24] E. B. Gledzer, "Systems of hydrodynamic type admitting two quadratic integrals of 
motion", Sov. Phys. Dokl. 18, 216 (1973). 

[25] M. Yamada and K. Ohkitami, "Lyapunov spectrum of a chaotic model of three- 
dimensional turbulence", J. Phys. Soc. Jpn 56, 4210 (1987). 

[26] S. Grossmann and D. Lohse, "Intermittency exponents", Europhys. Lett. 21, 201 (1993). 

[27] J. Eggers and S. Grossmann, "Anomalous turbulent velocity scaling from the Navier- 
Stokes equation", Phys. Lett. A 156, 444 (1991). 

[28] M. Farge, "Wavelet transforms and their applications to turbulence" , Ann. Rev. Fluid 
Mech. 24, 395 (1992). 

[29] C. Meneveau, "Analysis of turbulence in the orthonormal wavelet representation", J. 
Fluid Mech. 232, 469 (1991). 

[30] P. L. Garrido, J. L. Lebowitz, C. Maes, and H. Spohn, "Long-range correlations for 
conservative dynamics", Physical Rev. A 42, 1954 (1990). 

[31] C. W. Gardiner, "Handbook of Stochastic Methods", (Springer, Berlin, 1983). 

[32] R. H. Kraichnan, "Helical turbulence and absolute equilibrium", J. Fluid Mech. 59, 745 
(1973). 

[33] S. Orszag, "Lectures on the Statistical Theory of Turbulence" in: Fluid Dynamics 
1973, Les Houches Summer School of Theoretical Physics, edited by R. Balian and 



39 



J. L. Peube, (Gordon and Breach, 1977). 

[34] H. Tennekes and J. L. Lumley, A First Course in Turbulence, (MIT Press, 6th ed., 
1980). 

[35] S. Grossmann and D. Lohse, "Intermittency in the Navier-Stokes dynamics", Z. Phys. 
B 89, 11 (1992). 

[36] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes, 
(Cambridge University Press, 2nd edition, 1986). 

[37] R. M. Kerr and E. D. Siggia, "Cascade model of fully developed turbulence", J. Stat. 
Phys. 19, 543 (1978). 

[38] R. H. Kraichnan, "On Kolmogorov's inertial-range theories", J. Fluid Mech. 62, 305 
(1974). 

[39] C. Uhlig and J. Eggers, "Singularities in a cascade model of the Euler equation", 
manuscript in preparation (1996). 

[40] C. Uhlig and J. Eggers, "A stochastic model of turbulence on a fractal tree" , unpublished 
manuscript (1994). 

[41] A. J. Majda, "Vorticity, turbulence, and acoustics in fluid flow", SIAM Review 33, 349 
(1990). 

[42] P. Olla, "A closure model for intermittency in three-dimensional incompressible turbu- 
lence", Phys. Fluids 7, 1598 (1995). 

[43] E. D. Siggia, "Origin of intermittency in fully developed turbulence", Phys. Rev. A 15, 
1730 (1977). 

[44] M. H. Jensen, G. Paladin, and A. Vulpiani, "Intermittency in a cascade model for 
three-dimensional turbulence", Phys. Rev. A 43, 798 (1991). 



40 



[45] R. Benzi et al, "Scaling and dissipation in the GOY shell model", Phys. Fluids 7, 617 
(1995). 

[46] R. Benzi, L. Biferale, and G. Parisi, "On intermittency in a cascade model for turbu- 
lence", Physica D 65, 163 (1993). 



41 



